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' Abstract. We discuss the Schrodinger equation in presence of quaternionic potentials. 

, The study is performed analytically as long as it proves possible, when not, we resort to 

numerical calculations. The results obtained could be useful to investigate an underlying 
^ l' quaternionic quantum dynamics in particle physics. Experimental tests and proposals to 

U ■ observe quaternionic quantum effects by neutron interferometry are briefly reviewed. 
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I. INTRODUCTION 

After the classical mathematical and physical works on foundations of quaternionic quantum mechan- 
ics 0-^, there has been, in recent years, a widely interest in formulating quantum theories by using 
the non commutative ring of quaternions |p[-p^ . Some of the main results coming out from the use of 
new algebraic structures in particle physics are reviewed in the books of Dixon and Giirsey [ p^ . 
For a detailed discussion of quaternionic quantum mechanics and field theory we quote the excellent 
book of Adler ||l|l. 

The present paper has grown from an attempt to understand the experimental proposals [l^ 
|l8| and theoretical discussions ||l9|-|2l|l underlying the quaternionic formulation of the Schrodinger 
equation. The main difficulty in obtaining quaternionic solutions of physical problem is due to the 
fact that, in general, the standard mathematical methods of resolution break down. In the last years, 
some of these problems have been overcome. In particular, the discussion of quaternionic eigenvalue 
equations ||2^ and differential operators ||23|] is now recognized quite satisfactory. On the other hand, 
physical interpretations of quaternionic solutions represent a more delicate question . In discussing 
the Schrodinger equation what is still lacking is to understand the role that quaternionic potentials 
could play in quantum mechanics and where deviations from the standard theory would appear. 

The earliest experimental proposals to test quaternionic deviations from complex quantum me- 
chanics were made by Peres [|l6| who suggested that the non commutativity of quaternionic phases 
could be observed in Bragg scattering by crystal made of three different atoms, in neutron interfer- 
ometry and in meson regeneration. In 1984, the neutron interferometric experiment was realized by 
Kaiser, George and Werner p7| . The neutron wave function traversing slabs of two dissimilar materi- 
als (titanium and aluminumjshould experience the non commutativity of the phase shifts when the 
order in which the barriers are traversed is reversed. The experimental result showed that the phase 
shifts commute to better than one part in 3 x 10'*. To explain this null result, Klein postulated psl 
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that quaternionic potentials act only for some of the fundamental forces and proposed an experiment 
for testing possible violations of the Schrodinger equation by permuting the order in which nuclear, 
magnetic and gravitational potentials act on neutrons in an interferometer. 

The first theoretical analysis of two quaternionic potential barriers was developed by Davies and 
McKellar pl| |. In their paper, by translating the quaternionic Schrodinger equation into a pair of 
coupled complex equations and solving the corresponding complex system by numerical methods, 
Davies and McKellar showed that, notwithstanding the presence of complex instead of quaternionic 
phases, the predictions of quaternionic quantum mechanics differ from those of the usual theory. In 
particular, they pointed out that differently from the complex quantum mechanics prediction, where 
the left and right transmission amplitudes, i^, and t^, are equal in magnitude and in phase, in the 
quaternionic quantum mechanics only the magnitudes and \ta\ are equal. So, the measurement of 
a phase shift should be an indicator of quaternionic effects and of space dependent phase potentials. 
However, this conclusion leads to the embarrassing question of why there was no phase change in 
the experiment proposed by Peres and realized by Kaiser, George and Werner. To reconcile the 
theoretical predictions with the experimental observations, Davies and McKellar reiterated the Klein 
conclusion and suggested to subject the neutron beam to different interactions in permuted order. In 
the final chapter of the Adler book we find and intriguing question. Do the Kayser and colleagues 
experiment, and the elaborations on it proposed by Klein actually test for residual quaternionic effects? 
According to the non relativistic quaternionic scattering theory developed by Adler |l5[| the answer is 
clearly no. Experiments to detect a phase shift are equivalent to detect time reversal violation, which 
so far has not been detectable in neutron-optical experiments. 

In this paper, after a brief introductory discussion of quaternionic anti-self-adjoint operators, 
stationary states and time reversal invariance, we study the phenomenology of quaternionic one- 
dimensional square potentials. The j-k part of these potentials is treated as a perturbation of the 
complex case. We show that there are many possibilities in looking for quaternionic deviations from 
the standard (complex) theory. Nevertheless, in particular cases, we have to contend with quaternionic 
effects which minimize the deviations from complex quantum mechanics. With this paper, we would 
like to close the debate on the role that quaternionic potentials could play in quantum mechanics, but 
more realistically, we simply contribute to the general discussion. 



II. QUATERNIONIC SCHRODINGER EQUATION 

In the standard formulation of non-relativistic quantum mechanics, the complex wave function (p{r, t), 
describing a particle without spin subjected to the influence of a real potential V{r,t), satisfies the 
Schrodinger equation 

dMr,t)^j- [|iv^~F(r,i)] ^{r,t) . (1) 
In quaternionic quantum mechanics p5[ |, the anti-self-adjoint operator 

can be generalized by introducing the complex potential W{r, t) = t)\ exp[z0(r, t)], 

The anti-hermiticity is required to guarantee the time conservation of transition probabilities. As 
a consequence of this generalization for the anti-self-adjoint Hamiltonian operator, the quaternionic 
wave function ^{r,t) satisfies the following equation 

dMr, t) = { i [ 1^ - V{r, t)]+j. W{r, t) } <I>(r, t) . (2) 

Exactly as in the case of the standard quantum mechanics, we can define a current density 



^=2^ [(V<i>)^$-<I>^V<i>] 
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and a probability density 



Due to the non commutativity nature of quaternions, the position of the imaginary unit i in the 
current density is fundamental to obtain the continuity equation 



d,p + V ■J = 



(3) 



A. STATIONARY STATES 

The quaternionic Schrodingcr equation in presence of time-independent potentials 

[V{r),\Wir)\,e{r)] 

reads 



i-^V^-V{r)]+j.W{r)}<i>ir,t) . 



(4) 



The quaternionic stationary state wave function 

on the condition that \l/(r) be solution of the time-independent Schrodinger 



is solution of Eq 
equation 



-iV{r)+jW{r) ^'(r) + *(r) i £; = 



(5) 



Eq. represents a right complex eigenvalue equation on the quaternionic field |22 

Al'"'{r)-^{r) ^-^{r)iE . 

The allowed energies arc determined by the right complex eigenvalues X = i E of the quaternionic 
linear anti-sclf-adjoint operator yip" (r). The stationary state wave functions are particular solutions 
of Eq. (^). More general solutions can be constructed by superposition of such particular solutions. 
Summing over various allowed values of E, we get 



$(r,t) = ^ «'(r) cxp[- j.Et]q^ 



(6) 



where Qe arc constant quaternionic coefficients. The summation may imply an integration if the energy 
spectrum of E is continuous. 



B. TIME REVERSAL INVARIANCE 

From Eq. (|), wc can immediately obtain the time-reversed Schrodingcr equation 



2m 



(7) 



In complex quantum mechanics the ^-conjugation yields a time-reversed version of the original 
Schrodingcr equation. In quaternionic quantum mechanics there does not exist a universal time re- 
versal operator | [T5[ |. Only a restricted cla.ss of time-independent quaternionic potentials, i.e. 

W{r) = \W{r)\ exp[ie] , 

is time reversal invariant. For these potential, 

^T:{r,-t) = u^{r,t)u , u^kcxp[i9] . (8) 

For complex wave functions, we recover the standard result ^Ti^, —t) = $*(r, t). 
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C. ONE-DIMENSIONAL SQUARE POTENTIALS 

In solving the quaternionic Schrodinger equation, a great mathematical simplification results from the 
assumption that the wave function and the potential energy depend only on the x-coordinatc. 

i ^ ^{x) + [j W{x) - i V{x) ] *(a;) + iE^Q . (9) 

We shall consider one-dimensional problems with a potential which is pieced together from a number 
of constant portions, i.e. square potentials. In the potential region 

[V- \w\,e] 

the solution of the second order differential equation (^) is given by |Q 

^{x) {exp [z^'"'^!"'!:;;] Ci + exp [ - a; ] cj + 

w^'i"'!^^ {exp[z^'^'i'*'ix] C3 + exp[-z^-'''i"'ia;] c,} , (10) 

where Ci^....4 are complex coefficients determined by the boundary conditions. 



and 



, _ , IM^I cxp[ie] \ E;\w\,e _ [ • i \W\ cxp[-if)] 



In the free potential region, the solution reduces to 

= exp [ i | 2^] Ci + cxp [ — « f 2^] 0^+ j {cxp [fa;] c., + exp [ — fa;] C4} , 



where p = \/2mE. For scattering problems with a wave function incident from the left on quaternionic 
potentials, we have 

(.t) = exp[i I a; ] + r exp[ — i | x] + j r exp[ fa;], (11) 

where |r|^ is the standard probability of reflection and |f exp[ f a;]|^ represents an additional evanes- 
cent reflection. 



III. TIME REVERSAL INVARIANT (TRI) POTENTIAL BARRIER 

Let us consider the TRI potential 

[V{x); \W{x)\,0] . 

In Eq. (|^), the space-independent phase 6 can be removed by taking the transformation 

^{x) -> exp [i f ] ^{x) exp [-i f] . (12) 

Under this transformation 

Reflection and transmission probabilities do not change (actually the exponential exp[— i 9] can be 
absorbed in the complex coefficients C3.4). So, without lost of generality, we can discuss the quaternionic 
Schrodinger equation in presence of the square potential 

[V{x);\W{x)\] 

which has the following shape 
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Region I_ 



Region II_ 



Region III 



Region 



Region I+ 



[0;0] 



\V: 01 



[V; W] 



[V;0] 



The particle is free for x < —a, where the solution is given by Eq. and for x > a, where the 

solution is 



'^+(x) — t exp[i I a;] + ji exp[— | a;] 



(13) 



In Eqs. (11) and (|13|), wc have respectively omitted the complex exponential solutions cxp[ — ^ x] and 



exp[ ^ X ] because they are in conflict with the boundary condition that ^(x) remain finite as a; — > — oo 
and X — > +00. In order to determine the complex amplitudes r, f and we match the wave function 
and its slope at the discontinuities of the potential (see appendix). 

By using the continuity equation, we can immediately obtain the standard relation between the 
transmission and reflection coefficients, t and r. In fact, Eq. (|^) implies that the current density 

^ = [5xW(x)] I *(x) - ^(z) I d^^{x) ] 

has the same value at all points x. In the free potential regions, the probability current densities are 



given by ^_ = f ( 1 - \r\^ ) and 



Consequently, we find 



(14) 



In Figs. 1 and 2, we plot the transmission probability |t|^ as a function of i?[eV] for quaternionic 
potentials of different widths and heights (we assume that the incident particle is an electron). The 
presence of a quaternionic perturbation potential modifies the shape of the (complex) transmission 
probability curve. We have a reduction of the transmission probability. The presence of an inflection 
point is evident by increasing the width and the height of quaternionic potentials. Fig. 2 shows the 
transmission probability for the quaternionic potential 

V + j\W\ = (2.0 + j 1.5) eV 

and the complex comparative barrier [fol 



Z = VV^ + lW= 2.5eV . 
The wave numbers for these potentials are 



E > ^V^ + \W\^ 



\W\<E < ^/W+Wf 



E < \W\ 



v,w 

z 



v.w 



Z± 



\W[- 



V 



2m 



1/2 



\W\^ - E 



2m 



E 



where ip = arctan |^\/ (\W\'^ — E^) /F^j . For small quaternionic perturbations, the complex comparative 
barrier Z represents a good approximation of the quaternionic potential + j In this case, 



V I 



1 m 

2 y2 



and 
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The anti-sclf-adjoint operators corresponding to the quaternionic potential F + j and to the 
complex barrier Z are respectively 



Al 



V 



2m 



'jWir,t) and = 
By using these operators, we can write down two complex wave equations 
[Ay^'Y = -E^- and [Alf ^(x) = -E' 

The complex operators 



(15) 



\Al 



V 



\W\^ 



V - \W\' 



and 



[Al 



2m 



\w\ 



can be now easily compared. The difference is due to the factor which multiplies V^. Thus, complex 
comparative barriers only represent a first approximation to quaternionic potentials. In general, we 
have to consider the pure quaternionic potential, j \ W\, as a perturbation effect on the complex barrier 
V. 

Deviations from (complex) quantum mechanics appear in proximity of the complex barrier V when 
a quaternionic perturbation is turned on. Actually, in quaternionic quantum mechanics wc find an 
additional evanescent probability of transmission, that is This probability as a function of E[eV] 
is drawn in Fig. 3 for different values of x. 

To conclude the discussion of quaternionic one-dimensional time invariant potentials, we analyze 
the transmission probability \t\'^ as a function of the width of complex and quaternionic potentials. In 
Fig. 4, we plot the transmission probability for critical values of E. For E > \/V^ + the mini- 

mum value of the transmission probability oscillation decreases when the quaternionic perturbation 
increases. In Fig. 5, we compare complex and pure quaternionic potentials covering the same area, 
aV = h Clear deviations from complex quantum mechanics appear. 

IV. TIME REVERSAL VIOLATING (TRV) POTENTIAL BARRIER 

Let us modify the previous potential barrier by introducing a time reversal violating space-dependent 
phase d(x). Wc shall consider, for the region III, the following cases: 



Region IIIo 


Region IIIo 


Region IIIo 


Region III^ 


\V- |M^|,0] 


\V- |W^|,0] 


\V- \W\,d\ 


\V- \W\,Q\ 


Region IIIo 


Region IIIo 


Region lllg 


Region IIIq 


\V- |PF|,0] 


\V- \W\,B\ 


{V- \w\,e\ 


\V- |M^|,0] 


Region IIIj, 


Region IIIj, 


Region IIIo 


Region IIIq 


[F; \W\,d\ 


\V- \W\,d\ 


{V- |I^|,0] 


\V- |PF|,0] 



T 







As remarked in the introduction, quaternionic deviations from complex quantum mechanics could be 
observed by considering left and right transmissions through the same quaternionic potential barrier. 
The left transmission [x < —a) for the quaternionic potential of height \W\ and phase 



Q{x) 




(16) 
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is obviously equivalent to the right transmission [x > a) for the quaternionic potential of height \W\ 
and phase 

(9 -b<x <0 
9{x) = I . (17) 

[ < x < 

By using the transformation (^2|), we can replace the phase ( p7| ) by 

[ - 6 < .T < 
6{x) = <^ (18) 
[ -61 < a; < 6 

Thus, the plot of the transmission coefficient as a function of 9\n\ is a valid indicator of possible 
deviations from complex quantum mechanics. Symmetric curves (around the point 9\k\ ~ 1) shall 
imply no difference between left and right transmission through the same quaternionic barrier. In 
Figs. 6, 7 and 8, we draw the transmission probability, and the absolute value of the transmission 
coefficient, |j4r(7(t)|, as a function of the phase 9\k\. Qualitative deviations for complex quantum 
mechanics appear for asymmetric time violating potentials. It is also interesting to note that by 
increasing the phase [9\t^\ — *■ 1), quaternionic perturbation effects are minimized. For the convenience 
of the reader we explicitly give (see Tables 1, 2 and 3) the transmission probability |t|^ and the 
transmission coefficient t for different values of the potential phase 9 and the electron energy E. 



VI. CONCLUSIONS 

Very little progress in mathematical understanding of quaternionic analysis and algebra have often 
created (and sometimes justified) a distrust feeling to quaternionic formulations of physical theories. 
From our point of view, the recent increasing improvement of the mathematical structures involved 
in the quaternionic quantum mechanics could result in a rapid progress in this subject. 

The usefulness of quaternions (and, more in general, Clifford algebras) to unify algebraic and geo- 
metric aspects in discussing special relativity, Maxwell and Dirac equations is universally recognized. 
Nevertheless, notwithstanding the substantial literature analyzing quaternionic physical theories, a 
strong motivation forcing the use of quaternions instead of complex numbers is lacking. The experi- 
mental proposals of Peres , the theoretical analysis of Davics and McKcUar (l^,^ and the detailed 
and systematic development of quaternionic quantum mechanics in the Adler's book [jisl surely rep- 
resent the milestone in looking for quaternionic deviations from complex quantum mechanics. 

In this paper, we have presented a complete phenomenology of the quaternionic potential bar- 
rier by discussing the time invariant and time violating case. Interesting features of quaternionic 
perturbation effects emerge in the transmission and reflection coefScients. The various graphs show 
how the quantum measurement theory may be affected by changing from complex to quaternionic 
systems. The present work represents a preliminary step towards a significant advance in understand- 
ing quaternionic potentials and in looking for their experimental evidence. An interesting discussion 
about quaternionic violations of the algebraic relationship between the six coherent cross sections of 
any three scatterers, taken singly and pairwise, is found in p^ . 

Quaternionic time violating potentials and quaternionic perturbations (wich minimize the devia- 
tions from complex quantum mechanics) could play an important role in the CP violating physics. A 
theoretical discussion based on the wave packet formalism will be necessary to analyse experimental 
tests based on kaon regeneration [|6[|2j] . For asymmetric potentials a non null signals of quaternionic 
(time violating) effects should be observed. We will try to develop the wave packet treatment in a 
later article. 
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APPENDIX. MACTHING CONDITIONS 
A. TRI POTENTIAL BARRIER 

The matching conditions for the TRI potential barrier imply 



/1~ 

r 

r 



= S[a,b;E;V,\W\] 



(19) 



where 



S [a, b; E; V, \W\] = A_ M^' Dl_^ [M^]"^ x 



S[III] 



and 
D_ 



S[II] 



diag { exp [i I a] , exp [ - 1 1 a] , exp [ | a] , exp [ | a] } 



A_ 

QlWl.S 



1 -i ■ 
1 i ■ 



1 
^ 



M 



1 [w^^i"'''"], 



diag { exp [ ; 

V ' 

1 

0*1 



1 1 

;v^,|w| s: 



1 
1/' 



I 1 r e:v,\w\ 1 r e:v,\w\ 1 r i ^ 

77J , exp " ' 77J , exp " ' 77 J , exp [ - z+' " ' r] \ j 



1 

0-f,' ' 



D+ = diag{ exp [i|a] , exp [i |a] , exp [~ fa] , exp [- |a] } . 

The complex limit is obtained by setting 6 = 0. In this case (S[III] = 1)5 [a, b; E; V, W] reduces to 

S [a; E; V] = D_ A_ AP' D^^^ [M'"]-^ A^ . 

By matrix algebra, we easily calculate the coefficients for reflection and transmission 

t = exp [ — 2 i I a ] { cosh [ 2 ''^ a ] + | %_ sinh [ 2 2^^'^ a ] } ^ , 
r = — i x+ sinh [2 a] t 
i = , 
f = , 



where = | z^^^ ± (| z^^^) 
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B. TRV POTENTIAL BARRIER 

The matrix S [a, b, c; E; V, \ W\, 9] is now expressed in terms of 

rs[iii„„,J :5[0,-&] X s[e,-b] 

S[III] = J S[III„,,J : S[0,c~b] X S[9,-2c] x S[0,c- 
is[III,,„J X S[0-b] 

where 
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Table 1 [Fig. 6] 

E^O.beV; V = 2\W\ = 2eY; a = 26 = 4c=lA 



[E; V, \W\] 


e 




t 


[0.5; 2,0] 


[0,0; 0,0] 


0.62596 


0.618647 -i 0.493189 


[0.5; 2, 1] 


[0, 0; 0, 0] 


0.612889 


0.605527 0.496212 


[0.5; 2, 1] 




■ - - 

, ' ' 6 ' 6 . 

- ■ - 

, ' 6 ' 6 ' . 

- --00 




0.613720 
0.613748 
0.613720 


0.606559-^0.495789 
0.606389 -i 0.496025 
0.606159-^0.496278 


[0.5; 2, 1] 




■ - - 

' ' 4 ' 4 

- ■ - 

' 4 ' 4 ' 

'- - ■ 

47 4 ' ' . 




0.614708 
0.614769 
0.614708 


0.60763-^0.495473 
0.607413 -i 0.495800 
0.607065 -i 0.496165 


[0.5; 2, 1] 




■ - - 

. ' " ' 3 ' 3 . 

- ■ - 

' 3 ' 3 ' 

- - ■ 

. 3 ' 3 ' " ' " . 




0.615996 
0.616100 
0.615996 


0.608983 -i 0.495113 
0.608749 -i 0.495504 
0.608291 -i 0.495962 


[0.5; 2, 1] 




■ - - 

■ i ■ IL ' 

'- - ■ 6 

. 2 ' 2 ' ' . 




0.619112 
0.619321 
0.619112 


0.612155-^0.494347 
0.611982 -i 0.494772 
0.611358 -i 0.495332 


[0.5; 2, 1] 


[0, 0; TT, -k] 
[0, tt; ^,0] 
[tt, tt; 0, 0] 


0.625369 
0.625792 
0.625369 


0.618021 -^0.493376 
0.618478-^0.493232 
0.618021 -^0.493376 
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Table 2 [Fig. 7] 

E = 1.5eY] V = 2\W\ = 2cY ; a = 2& = 4c=lA 



[E; V, \W\] 


e 




t 


[1.5; 2,0] 


[0,0;0,0] 


0.845474 


0.835936-^0.382994 


[1.5; 2, 1] 


[0,0;0,0] 


0.840310 


0.830647-^0.387733 


[1.5; 2, 1] 




0- - - 

- ■ 2L 

. ' 6 ' 6 ' . 

- --no 




0.840637 
0.840651 
0.840637 


0.831080-^0.387223 
0.830996-^0.387424 
0.830877 -i 0.387660 


[1.5; 2, 1] 




■ - - 
-■ - 

' 4 ' 4 ' 

- 21 ■ 

41 4 ' ' . 




0.841023 
0.841055 
0.841023 


0.831516-^0.386787 
0.831409 -i 0.387058 
0.831228 -i 0.387406 


[1.5; 2, 1] 




0- - - 

. ' " ' 3 ' 3 . 

- ■ 2: 

. ' 3 ' 3 ' " . 

- --00 

. 3 ' 3 ' " ' " . 




0.841526 
0.841581 
0.841526 


0.832061 -z 0.386266 
0.831948-^0.386580 
0.831708 -i 0.387024 


[1.5; 2, 1] 




0- - - 

■ 2: ■ 2: ' 

' 2 ' 2 ' 

2: 2: ■ 




0.842740 
0.842850 
0.842740 


0.833323 -i 0.385115 
0.833247 -i 0.385421 
0.832917 -i 0.385991 


[1.5; 2, 1] 


[0,0; TT, tt] 
[0, ^; TT, 0] 
[tt, tt; 0, 0] 


0.845158 
0.845378 
0.845158 


0.835583-^0.383353 
0.835836-^0.383087 
0.835583-10.383353 
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Table 3 [Fig.8] 

= 3 cV ; y = 2 \W\ = 2eV; a = 2fe = 4c=lA 



[E; V, \W\] 







t 


[3; 2, 0] 


[0, 0; 0, 0] 


0.925842 


0.915930 -i 0.294811 


[3; 2, 1] 


[0, 0; 0, 0] 


0.923710 


0.913674 0.298177 


[3; 2, 1] 




■ - - 

, ' ' 6 ' 6 . 

- ■ - 

, ' 6 ' 6 ' . 

- --00 




0.923843 
0.923850 
0.923843 


0.913869 -i 0.297802 
0.913822-^0.297959 
0.913757 -i 0.298147 


[3; 2, 1] 




n ■ - - 

- ■ - 

' 4 ' 4 ' 

- - ■ 

47 4 ' ' . 




0.924000 
0.924016 
0.924000 


0.914057-^0.297490 
0.913998-^0.297699 
0.913898-^0.297977 


[3; 2, 1] 




■ - - 

. ' " ' 3 ' 3 . 

- ■ - 

' 3 ' 3 ' 

- - ■ b 




0.924205 
0.924232 
0.924205 


0.914289 -i 0.297122 
0.914226 -i 0.297360 
0.914095 -i 0.297718 


[3;2, 1] 




■ - - 

- 

'- - ■ 6 

. 2 ' 2 ' ' . 




0.924699 
0.924753 
0.924699 


0.914820-^0.296317 
0.914776-^0.296542 
0.914596 -i 0.297006 


[3; 2, 1] 


[0, 0; TT, n] 
[0, tt; ^,0] 
[tt, tt; 0, 0] 


0.925681 
0.925789 
0.925681 


0.915736-^0.295142 
0.915873 -i 0.294900 
0.915736 -i 0.295142 
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